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ABSTRACT 

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence 
depends on whether magnetic fields can change their topology via reconnection in highly conducting 
media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence 
is reassuring that the magnetic field behavior in computer simulations and turbulent astrophysical 
environments is similar, as far as magnetic reconnection is concerned. This makes it meaningful to 
perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field 
in astrophysical environments. Our studies of magnetic field diffusion in turbulent medium reveal 
interesting new phenomena. First of all, our three-dimensional MHD simulations initiated with anti- 
correlating magnetic field and gaseous density exhibit at later times a de-correlation of the magnetic 
field and density, which corresponds well to the observations of the interstellar media. While earlier 
studies stressed the role of either ambipolar diffusion or time-dependent turbulent fluctuations for 
de-correlating magnetic field and density, we get the effect of permanent de-correlation with one 
fluid code, i.e. without invoking ambipolar diffusion. In addition, in the presence of gravity and 
turbulence, our three-dimensional simulations show the decrease of the magnetic flux-to-mass ratio as 
the gaseous density at the center of the gravitational potential increases. We observe this effect both 
in the situations when we start with equilibrium distributions of gas and magnetic field and when we 
follow the evolution of collapsing dynamically unstable configurations. Thus the process of turbulent 
magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores 
and violently collapsing supercritical entities. The increase of the gravitational potential as well as 
the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated 
final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes 
the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is 
expected to play an important role in star formation, from its initial stages of concentrating interstellar 
gas to the final stages of the accretion to the forming protostar. In addition, we benchmark our codes 
by studying the heat transfer in magnetized compressible fluids and confirm the high rates of turbulent 
advection of heat obtained in an earlier study. 

Subject headings: diffusion — ISM: magnetic fields — magnetohydrodynamics (MHD) — star: forma- 
tion — turbulence 



1. INTRODUCTION 

Astrophysical flows are known to be turbulent and 
magnetized. The specific role played by MHD turbulence 
in different branches of astrophysics is still highly de- 
bated, but it is generally regarded as important. In par- 
ticular, for the interstellar medium (ISM) and star forma- 
tion, the ro le of turbulence has been discussed in many 
revie ws (see[Elmegrecn fc Scaldl200ilMcKee fc Ostrikerl 
|2007| ). The opinion on the role of magnetic field in 
these environments vary from magnetic field being re- 
garded as absolutely dominant in the proc esses (see 
iTassis fc Mouschoviasl [20051: iGalli et all [20061) to mod- 
erately import ant, as in super- Al fvenic models of star 
formation (see iPadoanet all[200l . 

The vital question that frequently permeates these de- 
bates is the diffusion of the magnetic field in astrophys- 
ical fluids. The conductivity of most of the astrophys- 
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ical fluids is high enough to make the Ohmic diffusion 
negligible on the scales involved, which means that the 
"frozen-in" approximation is a good one for many as- 
trophysical environments. However, without considering 
diffusive mechanisms that can violate the flux freezing, 
one faces problems attempting to explain many observa- 
tional facts. For example, simple estimates show that if 
all the magnetic flux is brought together with the ma- 
terial that collapses to form a star in molecular clouds, 
then the magnetic field in a proto-star should be sev- 
eral orders of magnitude higher than the one observed in 
T-Tauri stars ( this is the "magnetic flux problem", see 
Galli ct al. 2006] and references therein, for example) . 

To address the problem of the magnetic field diffu- 
sion both in the partially ionized ISM and in molecu- 
lar clouds, researcher s usually appeal to the ambipola r 
diffusion concept (see lMestel fc Spitzed[T95llShulll983f ). 
The idea of the ambipolar diffusion is very simple and 
may be easily exemplified in the case of gas collapsing 
to form a protostar. As the magnetic field is acting on 
charged particles only, it does not directly affect neu- 
trals. Neutrals move under the gravitational pull but 
are scattered by collisions with ions and charged dust 
grains which are coupled with the magnetic field. The 
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resulting flow dominated by the neutrals will be unable to 
drag the magnetic field lines and these will diffuse away 
through the infalling matter. This process of ambipolar 
diffusion becomes faster as the ionization ratio decreases 
and therefore, becomes more important in poorly ionized 
cl oud cores. 

IShu et ail (|2006l ) have explored the accretion phase in 
low-mass star formation and concluded that there should 
exist an effective diffusivity about 4 orders of magnitude 
larger than Ohmic diffusivity in order to an efficient mag- 
netic flux transport to occur. They have argued that 
ambipolar diffusion could work, but only under rather 
special circumstances like, for instance, considering par- 
ticular dust grain sizes. In other words, at the moment it 
is unclear if ambipolar diffusion is really high enough to 
solve the magnetic flux transport problem in collapsing 
flows. 

Does magnetic field remain absolutely frozen-in within 
highly ionized astrophysical fluids? The answer to this 
question affects the description of numerous essential 
processes in the interstellar and intergalactic gas. 

Mag netic reconnection was appealed in ILazariar] 
( 2005) as a way of removing magnetic flux from gravitat- 
ing clouds, e.g. from star-forming clouds. That work re- 
ferred to th e reconnection model of lLazarian fc Vishniad 
(1999) and lLazarian et all ()2004l ) for the justification of 
the concept of fast magnetic reconnection in the presence 
of turbulence. T he advantage of the scheme proposed by 
lLazarian! (|2005h was that robust removal of magnetic flux 
can be accomplished both in partially and fully ionized 
plasma, with only marginal dependence on the ionization 
state of the ga s 4 . The concept of "reconnection diffusion" 
introduced in lLazar ian ( 20Q5|) is relevant to our under- 
standing of many basic astrophysical processes. In par- 
ticular, it suggests that the classical textbook description 
of molecular clouds supported both by hourglass mag- 
netic field and turbulence is not self-consistent. Indeed, 
turbulence is expected to induce "reconnection diffusion" 
which should enable fast magnetic field removal from the 
cloud. However, in the absence of numerical confirma- 
tion of the fast reconnection, the scheme of magnetic flux 
removal through "reconnection diffusion" as opposed to 
ambipolar diffusion stayed somewhat speculative. 

For tunately, it has b een recently shown numerically 
(see iKowal et al.l |2009|) that three-dimensional mag- 
netic reconnection in turbulent flu id follows the predic- 
tions of the reconnection model of lLazarian fe Vishniad 
(1999) and therefore, is fast. This naturally in- 
creased the interest to the "reconnection diffusion" (see 
lLazarian fc Vishniac Il200"9l ). Motivated by this fact, here 
we perform simulations aiming to gain understanding of 
the diffusion of ma gnetic field induced by turbulence. As 
IKowal et all (|2009D tested magnetic reconnection in fully 
ionized gas in the present paper we shall focus our efforts 
on one fluid MHD simulations. 

We shall compare "reconnection diffusion" with the 
ambipolar diffusion and discuss the effects of ambipolar 
diffusion qualitatively focusing on the comparison of our 
results on "reconnection diffusion" with t hose on "am- 
bipolar turbulent diffusion" described in Hcitsch et al. 

4 The rates were predicted to depend on the reconnection rate, 
which according to Lazarian et al. ( 20(3) very weakly depends on 
the ionization degree of the gas. 



(2004). The latter study reported the enhancement of 
ambipolar diffusion in the presence of turbulence, which 
raised the question of how important is the simultaneous 
action of turbulence and ambipolar diffusion and whether 
turbulence alone, i.e. without any effect from ambipolar 
diffusion, can equally well induce de-correlation of mag- 
netic field and density. 

Our work on "reconnection diffusion" should also be 
distinguished from the research on the de-correlation of 
magnetic field and density with i n com p ressib le turbulent 
fluctuations. ICho fc Lazarian! (2002, 2003) performed 
three-dimensional MHD simulations and reported the ex- 
istence of separate turbulent cascades of Alfven and fast 
modes in strongly driven turbulence as well as a cas- 
cade of slow modes driven by Alfvenic cascade. Slow 
modes in magnetically dominated plasma are associated 
with density perturbations with marginal perturbation 
of magnetic fields, while the same is true for fast modes 
in weakly magnetized or high beta plasmas. Naturally, 
these two modes de-correlate magnetic fields and density 
on the crossing time of the wave. This was the effect stud- 
ied in more detail in one-dimensional both analyticall y 
and numerically bv lPassot fc Vazquez-Semadenil (|2003f> . 
who stressed that the enhancements of magnetic field 
strength and density may correlate and anti-correlate in 
turbulent interstellar gas within the fluctuations and this 
can introduce the dispersion of the mass-to-flux ratios 
within the turbulent volume. Each of the fluctuations 
provide a transient change of the pointwise magnetiza- 
tion. In the absence of other effects, e.g. related to 
the thermal instability, the de-correlation is reversible. 
In comparison, the "reconnection diffusion" , similar to 
the ambipolar diffusion, deals with the permanent de- 
correlation of magnetic field and density. Acting alone, 
the "reconnection diffusion" increases entropy making 
magnetic field-density de-correlation irreversible. 

What are the laws that govern magnetic field diffu- 
sion in turbulent magnetized fluids? Could those affect 
our understanding of basic interstellar and star forma- 
tion processes? These are the questions that we address 
in this paper. 

In this study, we try to understand the diffusion of 
magnetic field in a couple of idealized models in the pres- 
ence of turbulence. We explore setups both with and 
without gravity and compare the diffusion of magnetic 
field with that of a passive scalar. In the context of star 
formation an important issue that we will address is an 
alternative way of decreasing the magnetic flux-to-mass 
ratio without appealing to ambipolar diffusion. We claim 
that since turbulence in astrophysics is really ubiquitous, 
our results should be widely applicable. We also perform 
simulations including turbulent heat diffusion, which al- 
low us to compare results obtained with our code with 
those in the literature. 

This paper is organized as follows. In Section 2, we 
draw the theoretical grounds about fast magnetic recon- 
nection. In Section 3, we describe the numerical code 
employed. In Section 4, we present the results concerning 
the diffusion of magnetic field in a setup without external 
gravitational forces. In Section 5, we present the results 
of our experiments of diffusion of magnetic field in the 
presence of a gravitational field. In Section 6, we discuss 
our results and compare with previous works. In Section 
7, we discuss the accomplishments and limitations of our 
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present study. In Section 8, we discuss our findings in the 
context of strong turbulence theory, and finally in Sec- 
tion 9, we summarize our conclusions. While our work is 
focused on the diffusion of magnetic fields, we address in 
the Appendix the heat transport in magnetized turbulent 
plasma. We confir m with higher resolution the results in 
iCho et aT] (|2003af) that the heat advection within turbu- 
lent flows in the presence of magnetic fields is very similar 
to that induced by hydrodynamic turbulence. 

2. FAST MAGNETIC RECONNECTION IN THE PRESENCE 
OF TURBULENCE AND MAGNETIC DIFFUSION DUE 
TO RECONNECTION 

The dynamical response of magnetic fields in turbu- 
lent fluids, as we discussed above, depends on the ability 
of magnetic fields to change their topology via recon- 
nection. We know from observations that magnetic field 
reconnection may be both fast and slow. Indeed, a slow 
phase of reconnection is necessary in order to explain the 
accumulation of free energy associated with the magnetic 
flux that precedes eruptive flares in magnetized coronae. 
Thus it is important to identify the conditions for the 
reconnection to be fast. Different mechanisms prescribe 
different necessary requirements for this to happen. 

The problem of magnetic reconnection is most fre- 
quently discussed in terms of solar flares. However, this is 
a general basic process underlying the dynamics of mag- 
netized fluids in general. If the magnetic field lines in a 
turbulent fluid do not easily reconnect, the properties of 
the fluid should be dominated by intersecting magnetic 
flux tubes which are unable to pass through each other. 
Such fluids cannot be simulated with the existing codes 
as magnetic flux tubes readily reconnect in the numerical 
simulations which are currently very diffusive compared 
to the actual astrophysical flows. 

The famous Sweet-Par ker model of reconnection 
(|Sweetl 119581: iParkerl 119581: see Figure [Q upper panel) 
produces reconnection rates which are smaller than the 
Alfven velocity by a square root of the Lundquist num- 
ber, i.e. by S 1 " 1 / 2 = (L^Va/^) -1 / 2 , where L x in this case 
is the length of the current sheet. Astrophysical values of 
S can be as large as 10 15 or 10 20 , thus this scheme pro- 
duces reconnection at a rate which is negligible for most 
of the astrophysical circumstances. If the Sweet-Parker 
were proven to be the only possible model of reconnec- 
tion, it would have been possible to show that MHD nu- 
merical simulations do not have anything to do with real 
astrophysical fluids. Fortunately, faster schemes of re- 
connection are available. 

The first mo del of fast reconnection proposed by 
iPetschekl (|1964[ ) assumed that magnetic fluxes get into 
contact not along the astrophysically large scales of 
L x , but instead over a scale comparable to the resis- 
tive thickness S, forming a distinct X-point, where mag- 
netic field lines of the interacting fluxes converge at a 
sharp point to the reconnection spot. The stability 
of such a reconnection geometry in astrophysical situ- 
ations is an open issue. At least for uniform resistiv- 
ity, this configuration was proven to be unstable and to 
revert to a Sweet-Parker configuration (Biskamp' 1 119861 : 
lUzdenskv fc Kulsrudll2000l) . 

Recent years have been marked by the progress in un- 
derstanding some of the key processes of reconnection 
in astrophysical plasmas. In particular, a substantial 
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Fig. 1. — Upper plot: Sweet-Parker model of reconnection. The 
outflow is limited by a thin slot A, which is determined by Ohmic 
diffusivity. The other scale is an astrophysical scale L x 2> A. Mid- 
dle plot: reconnection of weakly stochastic magnetic field according 
to LV99. The model that accounts for the stochasticity of magnetic 
field lines. The outflow is limited by the diffusion of magnetic field 
lines, which depends on field line stochasticity. Low plot: an in- 
dividual small-scale reconnection region. The reconnection over 
small patches of magnetic field determines the local reconnection 
rate. The global reconnection rate is subs tantially larger as man y 
independent patches come together. From Lazarian ct al. (2004). 



progress has been obtained by c onsidering reconnectio n 
in the presence of the Hall-effect (IShav et al.l [l998. 2004). 
The condition for which the Hall-MHD term becomes im- 
portant for the reconnection is that the ion skin depth 
<5i on becomes comparable with the Sweet-Parker diffusion 
scale <5sp. The ion skin depth is a microscopic charac- 
teristic and it can be viewed as the gyroradius of an ion 
moving at the Alfven speed, i.e. S lon — VaI^cu where 
uj C i is the cyclotron frequency of an ion. For the parame - 
ters of the ISM (see Table 1 in lDraine & Lazarianlll998l) . 
the reconnection is collisional (see further discussion in 
lYamada et all [20061 ) . 

To deal with both collisio nal and collisionless plasma 
lLazarian fc Vishniad (j 19991 . henceforth LV99) proposed 
a model of fast reconnection in the presence of weak tur- 
bulence when magnetic field back-reaction is extremely 
important. We stress that the LV99 model does not as- 
sume that the magnetic field can be easily bent by fluid 
motions, which is a usual assumption to "justify" the er- 
roneous concept of turbulent diffusivity in the mean field 
dynamo. 

The middle and bottom panels of Figure [1] illustrate 
the key components of LV99 model 5 . The reconnec- 

5 The cartoon in Figure [T] is an idealization of the reconnec- 
tion process as the actual reconnection region also includes re- 
connected open loops of magnetic field moving oppositely to each 
other. Nevertheless, the cartoon properly reflects the role of the 
three-dimensionality of the reconnection process, the importance 
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tion events happen on small scales An where magnetic 
field lines get into contact. As the number of indepen- 
dent reconnection events that take place simultaneously 
is 1^/ Am 3> 1 the resulting reconnection speed is not 
limited by the speed of individual events on the scale A|| . 
Instead, the constraint on the reconnection speed comes 
from the thickness of the outflow reconnection region A, 
which is determined by the magnetic field wandering in 
a turbulent fluid. The model is intrinsically three dimen- 
sional 6 as both field wandering and simultaneous entry 
of many independent field patches, as shown in Figure HJ 
are three-dimensional effects. In LV99 model the mag- 
netic reconnection speed becomes comparable with Va 
when the scale of magnetic field wandering A becomes 
comparable with L X J 

For a quantitative description of the reconnection, 
one should adop t a model of MHD turbulence (see 
Iroshnikovl 11963c iKraichnanl 119651: iDobrowohrv et al.1 
1980t IShebalin et all 119831: iMontgomerv fc Turned 119811: 
Higdonl [1984D . The most important for magnetic field 
wandering is the A l fvenic component. Adopting the 
IGoldreich fc Sridharl (fl995l henceforth GS95) scaling of 
the Alfvenic component of MHD turbulence extended to 
include the case of weak turbulence, LV99 predicted that 
the reconnection speed in a weakly turbulent magnetic 
field is 

Vr = VAil/L^^/VA) 2 (1) 

where the level of turbulence is parameterized by the 
injection velocity; the combination Va(vi/Va) 2 is the 
velocity of the largest strong turbulence eddies V^trong, 
i.e., the velocity at the scale at which the Alfvenic tur- 
bulence transfers from the weak to the strong regimes. 
Thus Equation ((T|) can also be rewritten as Vr = 
V^trongG/ ' Lx) 1 / 2 ; vi < Va and I is the turbulence injection 
scale. 

The scaling predictions given by Equation ((J) have 
been tested successfull y by three-dimensi onal MHD nu- 
merical simulations in iKowal et al.l (|2009f ). This stimu- 
lates us to adopt the LV99 model as a starting point for 
our discussion of magnetic reconnection. 

How can An be determined? In LV99 model, as many as 
LjVA_i_A|| localized reconnection events take place, each 
of which reconnects the flux at the rate V TCC ^ i OC ai/A_L, 
where V TCCi i oca i is the velocity of local reconnection 
events at the scale Aii. The individual reconnection 
events contribute to the global reconnection rate, which 
in three dimensions becomes a factor of L x /X^ larger, 
i.e., 

^rec, global ~ L x j\^ V^ C c. local- (2) 

The local reconnection speed, conservatively assuming 
that the local events are happening at the Sweet -Parker 
rate, can be easily obtained by identifying the local resis- 

of small-scale reconnection events, and the increase of the outflow 
region compared to the Sweet— Parker scheme. 

6 Bidimensional numerica l simulations of turbulent reconnection 
in lKulp a-Dvbcl et al. (2009) show that the reconnection is not fast 
in this case. 

7 Another process that is determined by magnetic field wander- 
ing is the diffusion of energetic particles perpendicular to the mean 
magnetic field. Indeed, the coefficient of diffusion perpendicular to 
the magnetic field in the Milky Way is just an order of unity less 
than the coefficient of d iffusion parallel to the magnetic field (see 
Giacalonc & Jokipii 1999, and rcfcrcncs therein). 



tive region #sp with X± and using the relations between 
Am and A_l that follow from the MHD turbulence model. 
The corresponding calculations in LV99 provided the lo- 
cal reconnection rate v/S 1-1 / 4 . Substituting this local re- 
connection rate in Equation ([2]) one estimates the global 
reconnection speed, if this speed were limited by Ohmic 
resistivity, which is larger than Va by a factor S 1 ^ 4 . As 
a result, one has to conclude that the reconnection does 
not depend on resistivity. 

The LV99 model is applicable to situations when 
plasma effects are included, e.g. the Hall-MHD effect, 
which increases effective resistivities for local reconnec- 
tion events. While the latter point is difficult to test di- 
rectly with existing plasma codes, e.g., with PIC codes, 
due to the necessity of simulating both plasma micro- 
physics effects as well as mac r ophys ical effects of mag- 
netic turbulence, IKowal et al.l (|2009f) simulated the ac- 
tion of plasma effects by parameterizing them via anoma- 
lous resistivities. The values of such resistivities are a 
steep function of the separation between oppositely di- 
rected magnetic field lines, which also determines the 
current separating magnetic fluxes. With anomalous re- 
sistivities, the structure of the fractal current sheet of the 
turbulent reconnection changed substantially, but no sig- 
nificant changes of the reconnection rate were reported, 
which agrees well with the theoretical expectations of 
LV99. Within the model the explanation of this stems 
from the fact that the reconnection is already fast (i.e., 
independent of resistivity) even when small-scale recon- 
nection events are mediated by the Ohmic resistivity, 
while the bottleneck for the reconnection process is pro- 
vided by magnetic field wandering. Thus the increase of 
the local reconnection rate does not increase the global 
reconnection speed. 

In this work, we address the problems which are rel- 
evant to the reconnection in a partially ionized, weakly 
turbulent gas. The c orresponding model o f reconnec- 
tion was proposed in lLazarian et al.l (|2004l henceforth 
LVC04). The extensive calculations summarized in Ta- 
ble 1 in LVC04 show that the reconnection for realistic 
circumstances varies from O.IVa to 0.03Va, i.e., is also 
fast, which should enable fast diffusion arising from tur- 
bulent motions. 

The fact that magnetic fields reconnect fast in turbu- 
lent fluids ensures that the large-scale dynamics that we 
can reproduce well with numerical codes is not compro- 
mised by the difference in reconnection processes in the 
computer and in astrophysical flows. This motivates our 
present study in which we investigate diffusion processes 
in turbulent magnetized fluids via three-dimensional sim- 
ulations. We feel that the issue of the dimension of the 
simulations is crucial for the turbulence dynamics and 
the reconnection and, as a result, for our final results. 

Visualization of turbulent diffusion of heat or any pas- 
sive scalar field is easy within the GS95 model, which can 
be interpreted as a model of Kolmogorov cascade perpen- 
dicular to the local direction of the magnetic field (LV99 
and more discussion in Section 8). The corresponding 
eddies are expected to advect heat similarly to the case 
of the hydrodynamic heat advection (see Appendix A). 
The corresponding visualization of the magnetic field dif- 
fusion is more involved. Every time that magnetic field 
lines intersect each other, they change their configura- 
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tion draining free energy from the system. In the pres- 
ence of self-gravity this may mean the escape of magnetic 
field, which is a "light fluid" from the self-gravitating 
gaseous "heavy fluid" . Naturally, if the turbulence gets 
very strong the system gets unbound and then the mixing 
of magnetic field and gas, rather than their segregation 
is expected. In what follows we test these ideas numeri- 
cally. 

3. MHD EQUATIONS AND THE NUMERICAL CODE 

The systems studied numerically in this work is de- 
scribed by the resistive MHD equations, assuming an 
isothermal equation of state: 



dp 

dt 



V • (pu) = 



(3) 



P \ — +u- V ) u= -c^V/7+(VxB)xB-pV* + f (4) 



- = Vx(uxB) 



?7oh! 



,V 2 B 



(5) 



plus the divergenceless condition for the magnetic field 
V • B = 0. The spatial coordinates are given in units of a 
typical length L*. The density p is normalized by a ref- 
erence density p„, and the velocity field u by a reference 
velocity f7*. The constant sound speed c s is also given 
in units of £7», and the magnetic field B is measured in 
units of U*y/4irp*. The uniform Ohmic resistivity 77011m 
is given in units of £/*L*. In our numerical calculations 
we will use both non-zero and zero values of ?7ohm- In the 
latter case, the calculations will include only numerical 
resistivity. Time t is measured in units of L^/U*. The 
external gravitational potential VE' is given in units of U%. 
The source term f is a random force term responsible for 
injection of turbulence. 

The above equations are solved inside a three- 
dimensional box with periodic boundary conditions. 
We use a shock-capturing Godunov -type scheme with 
an HLL solver (see, for exa mple, iKowal et al.l 120071 
iFalceta-Goncalves et al.l [20081 ). Time integration is per- 
formed with the Runge-Kutta method of second or- 
der. Unless we say explicitly the opposite, we assume 
?7ohm = 0. 

We employ an isotropic, non-helical, solenoidal, delta 
correlated in time forcing f . This forcing acts in a thin 
shell around the wave number kt = 2.5(2-7r/L), that is, 
the scale of turbulence injection l ln j that is about 2.5 
times smaller than the box size L (in all the simula- 
tions, we choose L = 1 in code units). In most of the 
experiments, the rms velocity Kms induced by turbu- 
lence in the box is close to unity (in code units). There- 
fore, in these cases, the turnover time of the energy- 
carrying eddies (or the turbulent timescale tturb) is 
iturb ~ ^inj/wturb ~ (L/2.5)/V lms 0.4 units of time 
in code units. 

We note that for our present purposes we use an one- 
fluid approximation, which does not include ambipolar 
diffusion. This choice is appropriate for approximating a 
fully ionized gas. One may argue that the code describes 
also the dynamics of partially ionized gas, but on the 
scales where ions and neutrals are strongly coupled, i.e,. 
on scales lar ger than the scale of a mbipolar diffusion (see 
discussion in lLazarian et al.1l2004f) . 



4. TURBULENT MAGNETIC FIELD DIFFUSION IN THE 
ABSENCE OF GRAVITY 

Observati ons of different regions of the diffuse ISM 
compiled by iTroland fc Heilesl |l986) indicate that mag- 
netic fields and density are not straightf orwardly cor- 
relate d. These observations motivated iHeitsch et al.l 
(2004) to perform 2.5-dimensional numerical calculations 
in the presence of b oth ambipolar diffus ion and turbu- 
lence. The results in IHeitsch et al.l ()2004f ) also indicated 
de-correlation of magnetic field and density 8 and one 
may wonder whether ambipolar diffusion is always re- 
quired to de-correlate magnetic field and density or if, 
otherwise, to what extent the co ncept of "turb u lent a m- 
bipolar diffusio n" introduced in IHeitsch et al.l (|2004l ) is 
useful (see also iZw eibcl 2002|) . To address these issues 
we performed three-dimensional simulations of magnetic 
diffusion in the absence of ambipolar diffusion effects. 

4.1. Numerical Setup 

The magnetic field is assumed to have initially only the 
component in the z-direction. The initial configuration 
of the magnetic and density fields are: 



B z (x,y) = B 



,'27T 

B\ cos ( —x 



2tt 

T 1 



(6) 



p(x, y)=p -ji {BoBi cos (?fx) cos (?fy) 

+0.5 [Si cos (^z) cos (^y)] 2 }, (7) 

where (x, y) = (0, 0) is the center of the x, y-plane. 
Boundary conditions are periodic. 

This initial magnetic field configuration has an uniform 
component Bq plus an harmonic perturbation of ampli- 
tude B\. The density field is distributed in such a way 
that the gas pressure, given by the isothermal equation 
of state p = c 2 s p, equilibrates exactly the magnetic pres- 
sure, giving a magneto-hydrostatic solution. We choose 
the parameters B\ = 0.3, po = 1 and c s = 1 in all our 
simulations. Figure [2] illustrates these initial fields when 
B = 1.0. 

We should remar k that the simulations presented in 
IHeitsch et al.l (|2004D do not start at the equilibrium, like 
ours. There is no pressure term in their equation for the 
evolution of the momentum of the ions to counterbalance 
the magnetic pressure. In addition, in their work the ion- 
density field is kept constant in time and space. 

Another differenc e between our setup and that in 
IHeitsch et al.l (|2004D is that our parameter B\ is assumed 
to be the same for all the models studied, and not a frac- 
tion of Bq. Also, the amplitude of the perturbation of 
the homogen eous component of th e density field is a free 
parameter in IHeitsch et all ([20041 ) . while here it is con- 
strained by the imposed equilibrium between gas and 
magnetic pressures. 

8 We feel that th e constrained geometry of the simulations in 
IHeitsch ct al. (2004) (where the magnetic field was assumed to be 
perpendicular to the plane of the two-dimensional turbulence, so 
that there were no reconnection) weakened the comparison of their 
set up with effects in the magnetized ISM, but this point is beyond 
the scope of our present discussion. 
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Fig. 2. — (x, j/)-plane showing the initial configuration of the z component of the magnetic field B z (left) and the density distribution 
(right) for the model B2 (sec Tablc[l}. The centers of the plots correspond to (x,y) = (0,0). 



In addition, we introduce a passive scalar field $ ini- 
tially identical to B z . The parameters of our relevant 
simulations are presented in Table [T] 

In these simulations we keep the random velocity ap- 
proximately constant, i.e., V Yms f=s 0.8 for all the mod- 
els after one time step. Therefore, all these models are 
slightly subsonic. 

4.2. Notation 

Hereafter, the quantities within brackets with sub- 
script "i? = 0.25L": (-)fl=o.25L> or simply "0.25": (-) 25 
will denote averages inside a cylinder with main-axis in 
the z-direction centered in the computational box, with 
radius R = 0.25L, while a subscript "z": (•) will denote 
an average over the z-direction. An overbar means the 
average of some quantity inside the entire box. 

4.3. Results 

Figure [3] shows the evolution of the amplitude of the 
mode that is identical to the initial harmonic perturba- 
tion of the magnetic field (i.e., the rms of the amplitude 
of the Fourier modes (k x ,k y ) — (±1, ±1)), for (B z ) z and 
(B z ) z I (p) z . Most right plot in Figure UM (Appendix) 
shows the evolution of the amplitude of the same mode 
for ($) and ($) / (p) z - All the curves presented were 
smoothed in order to make the visualization clearer. We 
see that the decay of the magnetic field occurs at a sim- 
ilar rate to that of the passive field. The mode decays 
nearly exponentially at roughly the same rate for most 
of the models. Only the model Bl (Bq = 0.5) exhibits a 
higher decay rate. This may be due to the large scale field 
reversals that are common in super- Alfvenic turbulence. 
Table [T] shows the fitted values (and the uncertainty) of 
?7t ur b in the curves corresponding to the evolution of the 
amplitude of the modes for (B z ) . The fitted curve is 
exp { — fc 2 ^turb^}, where k 2 = k 2 + k 2 = 2 is the square 
of the module of the corresponding wave-vector. We ob- 
serve that the decay of the amplitude of the modes of 
(B z ) z is not continuous but saturates at a value that is 
naturally maintained by the turbulence (in Figure 3 it 
occurs after about t = 6). 

The diffusion of B/p on l arge scales was also ob- 
served in |Hcitsch et al. (2004]) for two-fluid simulations 
and there it was associated with the difference between 
the velocity field of the ions and neutrals, at small scales. 
However, here we observe a similar effect, but in one-fluid 
simulations, which is suggestive that turbulence rather 
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Fig. 3. — Evolution of the rms amplitude of the Fourier modes 
(k x ,k y ) = (±1,±1) of (B z ) z (upper curves) and (B z ) z /(p) z 
(lower curves). The curves for (B z ) z were multiplied by a fac- 
tor of 10. All the curves were smoothed to make the visualization 
clearer. 



than the details of the microphysics are responsible for 
the diffusion. 

Left and center panels of Figure 2] shows the distribu- 
tion of (p) z versus (B z ) z for the model B2 (B = I .0) at 
the initial configuration (t — 0) and after 10 time steps. 
We see in this projected view that the initial magnetic 
flux-to-mass relation is quickly spread and, in contrast 
with the $ p distribution (see Figure [I5j in the Ap- 
pendix), we do not see any tendency for the magnetic 
field and density to become correlated. 

To give a quantitative measure of the evolution of 
the flux-to-mass relation in the models, let us consider 
(SB, Sp) , the correlation between fluctuations of the mag- 
netic field SB and density Sp, defined by 

yj J (B - B) 2 d 3 Xx Jf(p - p) 2 d 3 x 

Right panel of Figure S] shows the evolution of (SB,5p) 
(see right side of Figure [15] for the evolution of (60>,8p) 
which is similarly defined). Differently from the passive 
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TABLE 1 

Parameters of the Simulations in the Study of Turbulent Diffusion of Magnetic Flux without Gravity. 



Model 


Bo 




^rras 


^turb 


''hurb _ B z 


%urb - * 


Resolution 


Bl 


0.5 


0.4 


0.8 


0.5 


0.18(0.03) 


0.19(0.02) 


256 3 


B2 


1.0 


0.4 


0.8 


0.5 


0.09(0.02) 


0.14(0.02) 


256 3 


B3 


1.5 


0.4 


0.8 


0.5 


0.10(0.02) 


0.08(0.02) 


256 3 


B4 


2.0 


0.4 


0.8 


0.5 


0.13(0.02) 


0.11(0.02) 


256 3 


B21 


1.0 


0.4 


0.8 


0.5 


0.15(0.02) 


0.11(0.02) 


128 3 


B2h 


1.0 


0.4 


0.8 


0.5 


0.10(0.01) 


0.13(0.02) 


512 3 



CD 



1.6 
1.4 
1.2 

r i.o 

0.8 
0.6 
0.4 
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Fig. 4. — Distribution of (p) vs. (B z ) for model B2 (see Table [TJ, at t = (left) and t = 10 (center). Bight: correlation between 
fluctuations of the strength of the magnetic field (SB) and density (Sp). 



scalar field that quickly becomes correlated to the density 
field, the magnetic field keeps a residual anti-correlation 
with it. 

A more careful analysis of our results indicates that the 
correlation between magnetic field intensity and density 
depends on the Mach number M s . For example, when 
we calculate the correlation (SB, Sp) using the turbulent 
models of Table [5] (study of diffusion of passive scalar 
fields, see Appendix A), we find weak positive correla- 
tions for the supersonic models and negative correlations 
for the subsonic ones. These correlations increase with 
M s . Thus, the anti-correlation detected in Figure 2] can 
be due to the slightly subsonic regime of the turbulence. 
These correlations and anti-correlations at this level can- 
not be excluded by the observa tional data as discussed, 
e.g., in lTroland fc Heilesl (|1986l ). We shall address this 
issue in more detail elsewhere. 

To summarize, the results of Figures |3] and 0] suggest 
that the turbulence can substantially change the flux-to- 
mass ratio B / p without any effect of ambipolar diffusion. 
The diffusion of the magnetic flux occurs in a rate similar 
to the rate of the turbulent diffusion of heat (passive 
scalar), even for sub-Alfvenic turbulence. 

As remarked in Section 2, we should emphasize that 
the efficient turbulent diffusion of magnetic field that 
we are observing in the simulations above is due to fast 
magnetic reconnection because otherwise, if the tangled 
magnetic lines by turbulence were not reconnecting, then 
they would be behaving like a Jello-type substance and 
this would make the diffusive transport of magnetic flux 
very inefficient (contrary to what is observed in the simu- 
la tions) . The issue of m agnetic reconnection was avoided 
in lHeitsch et al.1 (|2004l ) due to the settings in which mag- 
netic field was assumed perpendicular to the plane of the 



fluid motions. Magnetic reconnection, however, is an ef- 
fect that is present within any realistic three-dimensional 
setup. 

4.4. Effects of Resolution on the Results 

To convince the reader that the above results are not 
being affected by numerical effects, we ran one of the 
models (model B2) with increased and decreased reso- 
lutions (models B2h and B21, respectively, see Table [1]). 
Figure [5] compares the same quantities presented in Fig- 
ure [3] for these models. We do not observe significant 
difference between them. Thus, we can expect that the 
results presented for the models with resolution of 256 3 
are robust. 

5. "RECONNECTION DIFFUSION" IN THE PRESENCE OF 
GRAVITY 

The "reconnection diffusion" of magnetic field in the 
absence of gravity can represent the magnetic field dy- 
namics in diffuse interstellar gas where cloud self-gravity 
is not important. In molecular clouds, clumps and accre- 
tion disks, the diffusion of magnetic field should happen 
in the presence of gravity. We study this process below. 

5.1. Numerical Approach 

In order to get an insight into the magnetic field diffu- 
sion in a turbulent fluid immersed in a gravitational po- 
tential, we have performed experiments in the presence 
of a gravitational potential with cylindric symmetry 
given in cylindrical coordinates (R, <f>, z) by: 
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*(R > i? max ) 
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Fig. 5. — Comparison between models of different resolution: B2, 
B21, and B2h (Table[T]l- It presents the same quantities as in Figure 
® 

where R = is the center of the (x, j/)-plane, and we 
fixed i?* = 0.1L and R max — 0.45L (i = 1 is the size of 
the computational box, as remarked in Section 3). We 
assume a relatively high value of in order to limit 
the values of the gravitational force and prevent the sys- 
tem to be initially Parker-Rayleigh-Taylor unstable. We 
assume an outer cut-off R ma ,x on the gravitational force 
to ensure the cylindrical symmetry while using periodic 
boundary conditions. 

In one class of experiments, we start the simula- 
tion with a magneto-hydrostatic equilibrium with j3 = 
-fgas/-Pmag = c 2 s p/(B 2 /8n) constant. The initial density 
and magnetic fields are, respectively, 

p(R) = po exp {(*(i? max ) - *(fl))/c2(l + /T !)} (11) 

B 2 (R)=c s ^/2(3- 1 p(R). (12) 

Figure [6] illustrates this initial configuration for one of 
the studied models (model C2, see Table [5]). 

We restricted our experiments to the trans-sonic case 
c s = 1 (in most of the experiments, we keep Vrms ~ 0.8, 
see Table [2]) . We also fixed po — 1. The turbulence is 
injected at t = and we follow the evolution of (B z )r 
and (p)r for eight time steps. Table [2] lists the parame- 
ters used for these experiments, p and B z represent the 
average of the density and magnetic field over the entire 
box. Va,i refers to the initial Alfven speed of the system. 
The rms velocity of the system V rms is measured after 
the turbulence is well-developed. 

We have also performed experiments starting out of 
equilibrium, with homogeneous fields: the system starts 
in free fall. We leave the system to evolve for eight time 
steps applying turbulence from the very beginning. For 
a comparison, we have also performed these experiments 
without turbulence. The initial uniform magnetic field 
is parallel to the z-direction for these models. Table [3] 
lists the parameters for these runs. The listed values of 
P refer to the initial conditions. 

Concerning the diffusion of the magnetic field, in or- 
der to provide a quantitative comparison between the 



models, we have also performed simulations with simi- 
lar initial conditions to the models presented in Tabled 
but without turbulence and with the explicit presence of 
Ohmic diffusivity ?7ohm in the induction equation. As 
these models have perfect symmetry in the z-direction, 
we simulated only a plane cutting the z-axis, that is, 
they are 2.5-dimensional simulations. We use a resolu- 
tion comparable to the turbulent three-dimensional mod- 
els. Table [4] lists the parameters for theses runs. We 
simulated a three-dimensional model equivalent to the 
model E7rl, and we found exact agreement in the time 
evolution of the magnetic flux distribution (not shown). 
Therefore, we can believe that in this case these 2.5- 
dimensional simulations give results which are equivalent 
to three-dimensional simulations. 

5.2. Results 
5.2.1. Evolution of the Equilibrium Distribution 

Top row of Figure [7] shows the evolution of (B Z ) Q 25 
(left), (p) 025 (middle), and (B z ) 25 / (p) 025 (right), nor- 
malized by the respective characteristic values inside the 
box (B z , p and B z /p), for the models CI, C2, and C3 
(j3 = 1). We compare the evolution of these quantities 
for different strengths of gravity A, maintaining the other 
parameters identical. The central magnetic flux reduces 
faster the higher the value of A. The flux-to-mass ratio 
has similar behavior. The other plots of Figure [7] show 
the profile of the quantities (B z ) z (upper panels), (p) 
(middle panels), and (B z ) z / (p) z (bottom panels) along 
the radius R, each column corresponding to a different 
value of A both for t = (in magneto- hydrostatic equilib- 
rium and constant /3) and for t = 8. We see the deepest 
decay of the magnetic flux toward the central region for 
the highest value of A at t = 8. 

In Figure [8j we compare the rate of the "reconnection 
diffusion" when we change the turbulent velocity and 
maintain the other parameters identical as in models C2, 
C4 and C5. An inspection of the left panel shows that 
the central magnetic flux (B z ) 25 decreases faster for the 
two highest turbulent velocities. Fluctuations are higher 
in the cases with higher velocity. The central density 
however, gets smaller for the highest turbulent velocity. 
This is explained by the fact that the dynamic pressure is 
higher for the largest velocities. The central flux-to-mass 
ratio (B Z ) Q 25 / (p) 25 decays for the two smallest veloc- 
ities. However, for the largest velocity, it is not clear if 
this ratio decreases or not. Looking at the middle graph 
of Figure [5] (bottom row), we see that the central den- 
sity decreases for the highest forcing. This is indicative 
that the turbulence driving overcomes the gravitational 
potential making the system less bound. 

Top row of Figure |9] compares the evolution of (B z ) 25 
(left), (p) 25 (middle), and (B z ) 025 / (p) 25 (right), nor- 
malized by the characteristic average values inside the 
box (B z , p and B z /p), for models C2, C6, and C7 with 
different /3. Both the central magnetic flux and the flux- 
to-mass ratio decreases faster for the less magnetized 
model (f3 — 3.3). The other plots of Figure |9] show 
the radial profile of the quantities (B z ) (upper panels), 
(p) z (middle panels), and (B z ) z / (p) z (bottom panels) 
for each model. We can again observe a lower value of 
the flux in the central region (relative to the external re- 
gions) for the highest values of (3 at the time step t = 8. 
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Fig. 6. — Model C2 (see Tabled- Top row: logarithm of the density field; bottom, row: B z component of the magnetic field. Left column: 
central xy, xz, and yz slices of the system projected on the respective walls of the cubic computational domain, in t = 0; middle and right 
columns: the same for t = 3 (middle) and t = 8 (right). 



TABLE 2 

Parameters for the Models with Gravity Starting at Magneto-hydrostatic Equilibrium with Initial Constant /3. 



Model p 1 Va 



CI 
C2 
C3 
CI 
C5 
C6 
C7 
C21 
C2h 



^turb 



??turb 



^turb/Krms'i 



nj 



Resolution 



1.0 
1.0 
1.0 
1.0 
1.0 
3.3 
0.3 
1.0 
1.0 



1.4 
1.4 
1.4 
1.4 
1.4 
0.8 
2.4 
1.4 
1.4 



0.6 
0.9 
1.2 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 



1.26 
1.52 
1.95 
1.52 
1.52 
2.40 
1.18 
1.52 
1.52 



1.59 
1.74 
1.98 
1.74 
1.74 
1.20 
2.66 
1.74 
1.74 



0.8 
0.8 
0.8 
1.4 
2.0 
0.8 
0.8 
0.8 
0.8 



0.5 
0.5 
0.5 
0.3 
0.2 
0.5 
0.5 
0.5 
0.5 



g 0.005 
« 0.01 
« 0.03 
i 0.10 - 0.20 
g 0.30 
S 0.02 
« 0.01 



§ 0.015 

w 0.03 

RJ 0.09 
0.18 - 0.36 

> 0.37 

S 0.06 

x 0.03 



256 J 
256 3 
256 3 
256 3 
256 3 
256 3 
256 3 
128 3 
512 3 



1 Initial parameter for the plasma: /? = Pgas/Pmag 

2 The parameter A for the strength of gravity (see Equations 9 and 10) is given in units of c 2 s L. 



The contrast between the central and the more exter- the three models, being higher for the more magnetized 
nal values for the flux-to-mass ratio is quite different for models. This is expected, as turbulence brings the sys- 
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Fig. 7. — Evolution of the equilibrium models for different gravitational potential. The top row shows the time evolution of (B z ) 25 /B z 
(left), {p) 25 Ip (middle), and ((B z ) 25 / (p) 25) / (p/ B z ) (right). The other plots show the radial profile of {B z ) z (upper panels), (p) z 
(middle panels), and {B z ) z / (p) (bottom panels) for the different values of A in t = (magneto-hydrostatic solution with /3 constant, see 
Tabled and t = 8. Error bars show the standard deviation. All models have initial = 1.0. 



tern in the state of minimal energy. The effect of varying 
magnetization in some sense is analogous to the effect 
of varying gravity. The equilibrium flux-to-mass ratio 
is larger in both the case of higher gravity and higher 
magnetization. The physics is simple, the lighter fluid 
(magnetic field) gets segregated from the heavier fluid 
(gas). 



All in all, we clearly see that turbulence substantially 
influences the quasi-static evolution of magnetized gas in 
the gravitational potential. The system in the presence 
of turbulence relaxes faster to its minimum potential en- 
ergy state. This explains the change of the flux-to-mass 
ratio, which for years was a problem to deal with invok- 
ing ambipolar diffusion. 
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Fig. 8. — Evolution of the equilibrium models for different turbulent driving. The top row shows the time evolution of (B z ) 2 $ I B z 
(left), (p)q 25 lp (middle), and ((B z ) 2 s I (p)q 25) / (PI Bz) (right). The bottom row shows the radial profile of {B z ) z (left), (p) z (middle), 
and {B z ) z / (p) (right) for each value of the turbulent velocity V r ms, in t = (magneto-hydrostatic solution with constant) and t = 8. 
Error bars show the standard deviation. All models have initial = 1.0. See Table [2] 



TABLE 3 

Parameters for the Models with Gravity Starting 
out-of-equilibrium, with initially homogeneous fields. 



Model 





V A ,i 


A 


P 


B z 




Resolution 


Dl 


1.0 


1.4 


0.9 


1.0 


1.41 


0.8 


256 3 


Dla 


1.0 


1.4 


0.9 


1.0 


1.41 


0.0 


256 3 


D2 


3.3 


0.8 


0.9 


1.0 


0.77 


0.8 


256 3 


D2a 


3.3 


0.8 


0.9 


1.0 


0.77 


0.0 


256 3 


D3 


0.3 


2.4 


0.9 


1.0 


2.45 


0.8 


256 3 


D3a 


0.3 


2.1 


0.9 


1.0 


2.45 


0.0 


256 3 


Dll 


1.0 


1.1 


0.9 


1.0 


1.41 


0.8 


128 3 


Dlh 


1.0 


1.4 


0.9 


1.0 


1.41 


0.8 


512 3 



5.2.2. Equilibrium Models: Comparison of Magnetic 
Diffusivity and Resistivity Effects 

In terms of the removal of the magnetic field from 
quasi-static clouds, does the effect of "reconnection dif- 
fusion" act similar to the effect of diffusion induced by 
resistivity? To address this question, we have performed 
a series of simulations with enhanced Ohmic resistivity 
(see models of Table 0} . 

In Figure [TU1 we compare the evolution of (B Z ) R (at 
different radius) for model C2 of Table [5] with similar 
resistive models without turbulence of Table 21 with dif- 
ferent values of Ohmic diffusivity 7/ohm- The decay seems 
initially faster and comparable with the highest value of 
VOhm (»;ohm = 0.05). But after this initial phase, the 



TABLE 4 

Parameters for the 2.5-dimensional Resistive Models with 
Gravity Starting with Magneto-hydrostatic Equilibrium 
and Constant /3. 



Model 





A 


VOhm 


Resolution 


ElrO 


1 





0.6 


0.005 


256 2 


Elrl 


1 





0.6 


0.01 


256 2 


E2rl 


1 





0.9 


0.01 


256 2 


E2r2 


1 





0.9 


0.02 


256 2 


E2r3 


1 





0.9 


0.03 


256 2 


E2r4 


1 





0.9 


0.05 


256 2 


E3rl 


1 





1.2 


0.01 


256 2 


E3r2 


1 





1.2 


0.02 


256 2 


E3r3 


1 





1.2 


0.03 


256 2 


E4r2 


1 





0.9 


0.10 


256 2 


E4r3 


1 





0.9 


0.20 


256 2 


E5r3 


1 





0.9 


0.30 


256 2 


E6rl 


3 


3 


0.9 


0.01 


256 2 


E6r2 


3 


3 


0.9 


0.02 


256 2 


E6r3 


3 


3 


0.9 


0.03 


256 2 


E7rl 





3 


0.9 


0.01 


256 2 



turbulent model (C2) seems to have a behavior similar 
to the resistive models with 77011m between 0.01 and 0.02. 

Figure \W\ (Appendix) compares the turbulent models 
CI, C3, C4, C5, C6, and C7 of Table M with similar 
resistive models of Table [4] After roughly one time step, 
the model CI (weaker gravitational field) seems to be 
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Fig. 9. — Evolution of the equilibrium models for different degrees of magnetization (plasma (3 = Pgas/Pmag). The top row shows the 
time evolution of {B z ) 25 /B z {left), (p) Q 25 /p (middle), and ({B Z ) Q 25 / (p) 25) /(P/B~z) (right). The other plots show the radial profile 
of (B z ) z (upper panels), (p) (middle panels), and (Bz) / (p) (bottom panels) for each value of 0, in t = (magneto-hydrostatic solution 
with j3 constant) and t = 8. Error bars show the standard deviation of the data. See Table f2] 



consistent with a value of i)ohm between 0.005 or lesser, 
while the model C3 (stronger gravitational field) seems to 
be consistent with ryohm ~ 0.03. These results show that 
the effective turbulent magnetic diffusivity is sensitive to 
the strength of the gravitational field. 

The resistive simulations with increasing ?7ohm values 
are more comparable with models with increasing tur- 



bulent velocity. The model C4 (with smaller turbulent 
velocity) seems to be consistent with the resistive model 
with ?7ohm = 0.10. The model C5 (with larger turbulent 
velocity) seems to be more comparable with the model 
with 77ohm ~ 0.30 (or higher), however, in this case we 
cannot associate a representative value of ?7ohm due to 
the very quick diffusion which occurs even before t = 1, 
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Fig. 10. — Comparison between the model C2 (turbulent diffusivity) and resistive models without turbulence (see Table[4]|. All the cases 
have analogous parameters. 



when the turbulence becomes well developed. 

The turbulent curve for the less magnetized model (C6) 
seems to follow the resistive curve with 77011m = 0.02, 
while the more magnetized model (C7) is comparable 
to the resistive model with ryohm ~ 0.01. This result 
indicates that the effective turbulent diffusivity is also 
sensitive to the strength of the magnetic field. 

In summary, the results above indicate a correspon- 
dence between the two different effects. In other words, 
the turbulent magnetic diffusion may mimic the effects of 
Ohmic diffusion of magnetic fields in gravitating clouds. 
However, we should keep in mind that the physics of tur- 
bulent diffusion and Ohmic resistivity is different. Thus 
this analogy should not be overstated. 

5.2.3. Evolution of N on- equilibrium Models 

Figure [Til shows the same set of comparisons as in Fig- 
ure M for the models Dl, D2, and D3 of Table[3]— these 
models have started out of the equilibrium with a homo- 
geneous density and magnetic field in a free fall system. 
Besides the runs with turbulence, we also present, for 
comparison, the evolution for the systems without tur- 
bulence (models Dla, D2a, and D3a). The strong os- 
cillations seen in the evolution of the central magnetic 
flux and density for these models (which are more pro- 
nounced in the models without turbulence) are acoustic 
oscillations, since the time for the virialization of these 
systems is larger than the simulated period. We note that 
the initial flux-to-mass ratio does not change in the cases 
without turbulence. We also observe similar trends as in 
Figure^ the higher the value of j3, the faster the decrease 
of the central magnetic flux relative to the mean flux into 
the box. We also note that the radial profile of the flux- 
to-mass ratio for the turbulent models crosses the mean 
value for the models without turbulence at nearly the 
same radius. This is due to the fact that the effective 
gravity potential in all these simulations acts up to this 
radius approximately. 

This set of simulations shows that the change of mass- 
to-flux ratio can happen at the time scale of the gravi- 
tational collapse of the system and therefore, turbulent 
diffusion of magnetic field is applicable also to dynamic 
situations, e.g., to the formation of supercritical cores. 



5.3. Effects of Resolution on the Results 

As in the case of the study presented in Section 4, we 
would like to know how the results shown in this section 
arc sensitive to changes in resolution. Again, we ran some 
models employing higher resolution and we inspected the 
changes in the results concerned. 

Figure Q2] compares the evolution of some of the quan- 
tities studied through this section for models C21 and 
C2h — which are identical to C2, except by the lower 
(C21) and higher resolution (C2h, see Table [2]). It shows 
no significant difference between these models. Figure [T2l 
also depicts the evolution of the same quantities for the 
models Dl, Dll, and Dlh. Both models have the same 
parameters as in model Dl, but model Dlh (Dll) has 
higher (lower) resolution (see Table [3]) . Again, we see no 
disagreement between the models. 

Therefore, the results presented in this section are not 
expected to change with an increase in resolution. 

5.4. Magnetic Field Expulsion Revealed 

Both in the case of equilibrium and non-equilibrium 
we observe a substantial change of the mass-to-flux ra- 
tio. Even our experiments with no turbulence injection 
confirm that this process arises from the action of turbu- 
lence. As a result, in all the cases with gravity the turbu- 
lence allows magnetic field to escape from the dense core 
which is being formed in the center of the gravitational 
potential. 

6. DISCUSSION OF THE RESULTS: RELATIONS TO 
EARLIER STUDIES 

Through this work we have perf ormed the compariso n 
of our results with the study by iHeitsch et al.l (|2004l) . 
Below, we provide yet another outlook of the connection 
of that stud y with the prese nt paper. We also discuss 
the work bv lShu et all (|2006l) . which was the initial mo- 
tivation of our study of the diffusion of magnetic field in 
the presence of gravity. 

6.1. Comparison with \Heitsch et al\ ^20 OA ): Ambipolar 
Diffusion Versus Turbulence and 2. 5- dimensional 
Versus Three-dimensional 

In view of the astrophysical impli cations, the compari - 
son between our results and those of lHeitsch et al.l (|2004[ ) 
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calls for the discussion on how ambipolar diffusion and 
turbulence intera ct to affect the magn etic field diffusiv- 
ity. In particular, Heit sch et al.l (|2004D claim th at a new 
proce ss "turbulent ambipolar diffusion" (see also lZweibell 
2002) acts to induce fast magnetic diffusivity. 

At the same time, our results do not seem to ex- 
hibit less magnetic diffusivity than those of Heitsch et al. 



( 2004) in spite of the fact that we do not have ambipolar 
diffusion. How can this be understood? We propose the 
following explanation. In the absence of ambipolar dif- 
fusion, the turbulence propagates to smaller scales mak- 
ing small-scale interactions possible. On the other hand, 
ambipolar diffusion affects the turbulence, increasing the 
damping scale. As a result, the ambipolar diffusion acts 
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in two ways, in one to increase the small-scale diffusivity 
of the magnetic field, in another is to decrease the tur- 
bulent small-scale diffusivity and these effects essentially 
compensate each other 9 . 

In other words, if we approximate the turbulent dif- 
fusivity by (l/3)V[nj£mj, where V5 n j and L; n j are the in- 
jection velocity and t he injection scal e for strong MHD 
turbulence (see LV99, lLazariar3l2006f) . respectively, the 
ambipolar diffusivity acting on small scales will not play 
any role and the diffusivity will be purely "turbulent". 
If, however, the ambipolar diffusion coefficient is larger 
than VinjLinj, then the Reynolds number of the steered 
flow may become small for strong MHD turbulence to 
exist and the diffusion is purely ambipolar in this case. 
We might speculate that this leaves little, if any, pa- 
rameter space for the "turbulent ambipolar diffusion" 
when turbulence and ambipolar diffusion synergetically 
enhance diffusivity, acting in unison. This point should 
be tested by three-dimensional two-fluid simulations ex- 
hibiting both ambipolar diffusion and turbulence. 

In view of our findings one may ask whether it is sur- 
prising to observe the "reconnection diffusion" of mag- 
netic field being independent of ambipolar diffusion. We 
can appeal to the fact well known in hydrodynamics, 
namely, that in a turbulent fluid the diffusion of a pas- 
sive contaminant does not depend on the microscopic dif- 
fusivity. In the case of high microscopic diffusivity, the 
turbulence provides mixing down to a scale l\ at which 
the microscopic diffusivity both, suppresses the cascade 
and ensures efficient diffusivity of the contaminant. In 
the case of low microscopic diffusivity, turbulent mixing 
happens down to a scale I2 "C h, which ensures that 
even low microscopic diffusivity is sufficient to provide 
efficient diffusion. In both cases the total effective diffu- 
sivity of the contaminant is turbulent, i.e. is given by the 
product of the turbulent injection scale and the turbu- 
lent velocity. This analogy is not directly applicable to 

9 A possible point of confusion is related to the difference of the 
physical scales involved. If one associates the scale of the reconnec- 
tion with the thickness of the Sweet— Parker layer, then, indeed, the 
ambipolar diffusion scale is much larger and therefore the recon- 
nection scale gets irrelevant. However, within the LV99 model of 
reconnection, the scale of reconnection is associated with the scale 
of magnetic field wandering. The corresponding scale depends on 
the turbulent velocity and is not small. 



ambipolar diffusion, as this is a special type of diffusion 
and magnetic fields are different from passive contami- 
nants. However, we believe that our results show that to 
some extent the concept of turbulent diffusion developed 
in hydrodynamics carries over (due to fast reconnection) 
to magnetized fluid. 

6.2. Transient De- correlation of Density and Magnetic 

Field 

Magneto-sonic waves are known to create transient 
changes of the density and magnetic field correlation. 
In the case of turbulence the situation is less clear, 
but the research in the field suggests that the de- 
composition of the turbulent motions into basic MHD 
modes is meanin gful and justified even for high am- 
pli tude motions dCho et al.l l2003af) . Thus the claim 
in iPassot fc Vazquez-Semadenil (|2003l) that even in the 
limit of ideal MHD, turbulence can transiently affect the 
magnetic field and density correlations is justified. How- 
ever, the process discussed in this paper is different in the 
sense that the de-correlation we describe here is perma- 
nent and it will not disappear if the turbulence dissipates. 
In a sense, as we showed above, "reconnection diffusion" 
is similar to the ambipolar and Ohmic diffusion. It is a 
dissipative diffusion process, which does require non-zero 
resistivity, although this resistivity can be infinitesimally 
small for the LV99 model of fast reconnection in the pres- 
ence of turbulence (see Section 2). 

6.3. Relation to \Shu et al\ i200(3l ): Fast Removal of 
Magnetic Flux During Star Formation 

As discussed in iShu et ail (|2006| ). the sufficiency of 
the ambipolar diffusion efficiency for explaining obser- 
vational data of accreting proto-stars is questionable. At 
the same time, they found that the required dissipation 
is about 4 orders of magnitude larger than the expected 
Ohmic dissipation. Thus they appealed to the hyper- 
resistivity concept in order to explain the higher dissipa- 
tion of magnetic field. 

We feel, however, that the hyper-res istivity idea is 
poorl y ju stified (see criticis m of it in Lazaria n et al.l 
120041 and iKowal et al.l 12009ft . At the same time, fast 
three-dimensional "reconnection diffusion" can provide 
the magnetic diffusivity that is adequate for fast remov- 
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ing of the magnetic flux. This is what, in fact, was 
demonstrated in the present set of numerical simulations. 

It is worth mentioning that, unlike the actual Ohmic 
diffusivity, "reconnection diffusion" does not transfer the 
magnetic energy directly into heat. The lion share of the 
energy is being released in the form of kinetic energy, 
driving turbulence (see LV99). If the system is initially 
laminar, this potentially result in flares of reconnection 
and the corresponding diffusivity. This is in agreement 
with LV99 scheme where a more intensive turbulence 
should induce more intensive turbulent energy injection 
and lead to the unstable feeding of the energy of the de- 
formed magnetic field. However, the discussion of this 
effect is bey ond the scope of the present paper. 

Similar to iShu et al.l (|2006l ). we expect to observe the 
heating of the media. Indeed, although we do not expect 
to have Ohmic heating, the kinetic energy released due to 
magnetic reconnection is dissipated locally and therefore 
we expect to observe heating in the medium. Our setup 
for gravity can be seen as a toy model representing the 
situation in IShu et all (|2006T) . In the broad sense, our 
work confirms that a process of magnetic field diffusion 
that does not rely on ambipolar diffusion is efficient. 

We accept that our setup assuming an axial gravita- 
tional field is a very simple and ignores complications 
that could arise from using a nearly spherical potential of 
the self-gravitating cloud. The periodic boundary condi- 
tions give super-stability to the system, and do not allow 
inflow (or outflow) of material/magnetic field as we ex- 
pect in a more realistic accretion process. However, our 
experiments can give us qualitative insights. They show 
that the turbulent diffusion of the magnetic field can re- 
move magnetic flux from the central region, leading to 
a lower fiux-to-mass ratio in regions of higher gravity 
compared with that of lower gravity. 

We chose parameters to the simulations such that the 
system is not initially unstable to the Parker-Rayleigh- 
Taylor (PRT) instability. Although the PRT instability 
could be present in real accretion systems and could help 
to remove magnetic field from the core of gravitational 
systems, its presence would make the interpretation of 
the results more difficult and we wanted to analyze only 
the turbulence role in the removal of magnetic flux. How- 
ever, it is possible that this instability had been also act- 
ing due to local changes of parameters due to the turbu- 
lent motion. To ensure that the transport of magnetic 
flux is being caused by injection of turbulence only, we 
stopped the injection after a few time steps in some ex- 
periments and left the system to evolve. When we did 
this, the changes in the profile of the magnetic field and 
the other quantities stopped. 

We showed that the higher the strength of the gravi- 
tational force, the lower the flux-to-mass ratio is in the 
central region (compared with the mean value in the com- 
putational domain). This could be understood in terms 
of the potential energy of the system. When the gravita- 
tional potential well is deeper, more energetically favor- 
able is the pile up of matter near the center of gravity, 
reducing the total potential energy of the system. When 
the turbulence is increased, there is an initial trend to 
remove more magnetic flux from the center (and conse- 
quently more inflow of matter into the center), but for 
the highest value of the turbulent velocity in our experi- 
ments, there is a trend to remove material (together with 



magnetic flux) from the center, reducing the role of the 
gravity, due to the fact that the gravitational energy be- 
came small compared to the kinetic energy of the system. 
Our results also showed that when the gas is less mag- 
netized (higher f3, or higher values of the Alfvenic Mach 
number Ma), the reconnection diffusion of magnetic flux 
is more effective, but the central flux-to-mass ratio rel- 
ative to external regions is smaller for more magnetized 
models (low /?), compared to less magnetized models. 
That is, the contrast B/p between the inner and outer 
radius is higher for lower (3 (or Ma)- 

If the turbulent diff usivity of magn etic field may ex- 
plain the results in IShu et ahl ([2006D . one may won- 
der whether one can remove magnetic field b y this way 
not on ly from the class of systems studied by IShu et all 
(2006), but also from less dense systems. For instance, it 
is frequently assumed that only ambipolar diffusion is im- 
portant for the evolution of su b-critical magnetized clouds 
(|Tassis fc Mouschoviasll2005l ). Our study indicates that 
this conclusion may be altered in the presence of tur- 
bulence. This point, however, requires further careful 
study, which is beyond the scope of the present paper. 
In the future, we intend to study a more realistic model, 
e.g., with open boundary conditions and more realistic 
gravitational potentials. 

7. TURBULENT MAGNETIC DIFFUSION AND 
TURBULENCE THEORY 

The concept of "reconnection diffusion" is related 
to the LV99 model of fast reconnection which makes 
use of the model of stro ng turbulence proposed by 
IGoldreich fc Sridharl (H995L henceforth GS95) the turbu- 
lence is being injected at the large scales with the in- 
jection velocity V \ n \ equal to the Alfven velocity Va (see 
iCho et aLll2003bl for a review) . The turbulent eddies mix 
up magnetic field mostly in the direction perpendicular 
to the local magnetic field thus forming a Kolmogorov- 
type picture in terms of perpendicular motions. Natu- 
rally, these eddies are as efficient as hydrodynamic ed- 
dies are expected in terms of heat advection. One also 
can visualize how such eddies can induce magnetic field 
diffusion. 

It is important to note that the GS95 model deals with 
motions with respect to the local rather than mean mag- 
netic field. Indeed, it is natural that the motions of the 
parcel of fluid are affected only by the magnetic field of 
the parcel and of the near vicinity, i.e., by local fields. 
At the same time, in the reference frame of the mean 
field, the local magnetic fields of different parcels vary 
substantially. Thus we do not expect to see a substantial 
anisotropy of the heat advection when Vjnj ~ Va- 

It was noted in LV99 that one can talk about turbu- 
lent eddies perpendicular to the magnetic field only if the 
magnetic field can reconnect fast. The rates of reconnec- 
tion predicted in LV99 ensured that the magnetic field 
changes topology over one eddy turnover period. If the 
reconnection were slow, the magnetic fields would form 
progressively complex structures consisting of unresolved 
knots, which would invalidate the GS95 model. The re- 
sponse of such a fluid to mechanical perturbations would 
be similar to "Jello" , making the turbulence-sponsored 
diffusion of magnetic field and heat impossible. 

What happens when Vj n j < Va? In this case the tur- 
bulence at large scales is weak and therefore magnetic 
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field mixing is reduced. Thus one may expect a partial 
suppression of magnetic diffusivity. However, as turbu- 
lence cascades the strength of interactions increases and 
at a scale L^(V^/Va) 2 th e turbulence gets strong. Ac- 
cording to iLazarianl (|2006| ) , the diffusivity in this regime 
decreases by the ratio of (V"i n j/V/i) 3 , with the eddies 
of strong turbulence playing a critical role in the pro- 
cess. When we compare the turbulent diffusivity Tfturb 
estimated for the sub-Alfvenic models described in Ta- 
ble [5] (see Section 5) with Li n jVt UT b(Vt ur b/VA) 3 , we find 
that th e values are rou ghly consistent with the predic- 
tions of lLazarianl (|2006f ) , although a more detailed study 
is required in t his regard. For instance, we know that 
ILazarianl (|2006t ) theory was not intended for high Mach 
number turbulence. 

All in all, we believe that the high diffusivity that 
we observe is related to the properties of strong mag- 
netic turbulence. While the latter is still a theory 
whic h is subject to intensive stud y (see jBoldvrevl I2005L 
200a iBeresnyak fc Lazarianll2006l t2009al|bl : iGogoberidzd 



2003), we believe that for the purpose of describing 



magnetic and heat diffusion the existing theory and the 
present model catch all the essential phenomena. 

8. ACCOMPLISHMENTS AND LIMITATIONS OF THE 
PRESENT STUDY 

8.1. Major Findings 

This paper presents several sets of simulations which 
deal with magnetic diffusion in turbulent fluids. Com- 
paring our result on magnetic diffusion and that of heat, 
we see many similarities in these two processes. Our nu- 
merical testing in the paper would not make sense if the 
astrophysical reconnection were slow. Indeed, the major 
criticism that can be directed to the work of turbulent 
diffusion of heat by iCho et all (|2003a| ) is that reconnec- 
tion in their numerical simulations was fast due to high 
numerical diffusivity. With the confi rmation of th e LV99 
model of turbulent reconnection by iKowal et al.l (2009) 
one may claim that astrophysical reconnection is also 
generally fast and the differences between the computer 
simulations and astrophysical flows are not so dramatic 
as far as the reconnection is concerned. 

The most important part of our study is the removal 
of magnetic fields from gravitationally bounded systems 
(see Section 5). Generally speaking, this is what one 
can expect on the energetic grounds. Magnetic field can 
be identified with a light fluid which is not affected by 
gravity, while the matter tends to fall into the gravita- 
tional potential 10 . Turbulence in the presence of mag- 
netic reconnection helps "shaking off" matter from mag- 
netic fields. In our simulations the gravitational energy 
was larger than the turbulent energy. In the case when 
the opposite is true, the system is expected to get un- 
bounded with turbulence mixing magnetic field in the 
same way it does in the absence of gravity (see Section 

It is important to note that in Section 5 we obtained 
the segregation of magnetic field and matter both in the 
case when we started with equilibrium distribution and 
in the case when the system was performing a free fall. In 
the case of non-equilibrium initial conditions the amount 

10 As a matter of fact, in our low simulations, which we did not 
include in the paper, we see clear signatures of the PRT instability. 



of flux removed from the forming dense core is substan- 
tially larger than in the case of the equilibrium magnetic 
field/density configurations (compare Figures l9l and [TTj) . 
Nevertheless, the flux removal happens fast, essentially 
in one turnover of the turbulent eddies. In comparison, 
the effect of numerical diffusion for the flux removal in 
our simulations is marginal, and this is testified by the 
constant flux-to-mass ratio obtained in the simulations 
without turbulence (see Figure [TTj) . 

What is the physical picture corresponding to our find- 
ings? In the absence of gravity turbulence mixes up 11 
flux tubes with different magnetic flux-to-mass ratios de- 
creasing the difference in this ratio. In the presence of 
gravity, however, it is energetically advantageous of flux 
tubes at the center of the gravitational potential to in- 
crease the mass-to-flux ratio. This process is enabled 
in highly conducting fluid by turbulence which induces 
"reconnection diffusion" . 

8.2. Applicability of the Results 

The diffusion of magnetic field in our numerical runs 
exhibits a few interesting features. First of all, according 
to Figure [4] one may expect to see a broad distribution 
of magnetic field intensity with density. This seems to 
be consistent with the m easurements of magnet ic field 
strength in diffuse media (|Troland fc Heileslll986l ) . 

The situation gets even more intriguing as we discuss 
magnetic field diffusion in the gravitational potential. It 
is tempting to appl y these results to star formation pro- 
cess (see studies bv lLeao et alJ l2009). There, molecular 
clouds are known to be either m agnetically supercrit i- 
cal or magnetically subcritical fsee lMestel fc R av 1985). 
If, however, the magnetic flux can be removed from the 
gravitating turbulent cloud in a timescale of about an 
eddy turnover time, then the difference between clouds 
with different initial magnetization becomes less impor- 
tant. The initially subcritical turbulent clouds can lose 
their magnetic flux via the turbulent diffusion to become 
supercritical. 

An important point of the turbulent diffusion of the 
magnetic field is that it does not require gas to be weakly 
ionized, which is the requirement of the action of the am- 
bipolar diffusion. Therefore, one may expect to observe 
gravitational collapse even of the highly ionized gas. 

8.3. Magnetic Field Reconnection and Different Stages 
of Star Formation 

"Reconnection diffusion" seems to be a fundamental 
process that accompanies all the stages of star forma- 
tion. Our work shows (Section 4) that three-dimensional 
diffusion of magnetic field provides a wide distribution 
of the mass-to-flux ratios with some of the fluctuations 
having this ratio rather high. We believe that the diffu- 
sion of magnetic field described here is one the reasons 
for creation of zones of s uper- Alfenic turbulenc e even for 
sub-Alfvenic driving (see iBurkhart et"alll2009| ). 

The regions of density concentration get gravitation- 
ally bound. One can associate such regions with GMCs. 

11 This mixing for Alfvenic modes happens mostly perpendicu- 
lar to the local magnetic field for sub-Alfvenic and trans-Alfvenic 
turbulence (LV99). For super- Alfvenic turbulence the mixing is 
essentially hydrodynamic at large scales and the picture with mo- 
tions perpendicular to the local magnetic field directio n is restored 
at small scales (see discussion in Lazarian ct al. 2004). 
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These entities are known to be highly turbulent and tur- 
bulent diffusion will proceed within them, providing a 
hierarchy of self-gravitating zones with different density 
and different mass-to-flux ratios. Some of those zones 
may be subcritical in terms of magnetic field and some 
of them may be supercritical. In subcritical magnetic 
cores the turbulent diffusion may proceed quasi-statically 
as we described in Section 5.2.1. and in the supercritical 
cores the turbulent diffusion may proceed as we described 
in Section 5.2.3. In both cases, we expect the removal 
of magnetic field from the self-gravitating cores. This 
process proceeds all the time, including the stage of the 
accretion disks. 

Indeed, in the turbulent scenario of star formation it 
usually is assumed that at the initial stages of star for- 
mation the concentration of mater ial happens due to ga s 
moving along magnetic field lines (jMestel fc P aris 1984). 
This one-dimensional process requires rather long times 
of the accumulation of material. In contrast, the "re- 
connection diffusion" allows for the much faster three- 
dimensional accumulation. 

What is the relative role of the ambipolar diffusion and 
the "reconnection diffusion" . This issu e requires f urther 
studies. It is clear from the study by IShu et aT] (|2006T ) 
that in some situations the ambipolar diffusion may be 
not fast enough to explain the removal of magnetic fields 
from accretion disks. This is the case when we claim that 
the "reconnection diffusion" should dominate. At the 
same time, in cores with low turbulence, the ambipolar 
diffusion may dominate the reconnection diffusion. The 
exact range of the parameters for one or the other process 
to dominate should be defined by future research. 

8.4. Unsolved Problems and Future Studies 

Our paper has a clearly exploratory character. For in- 
stance, to simplify the interpretation of our results we 
studied the concentration of material in the given grav- 
itational potential, ignoring self-gravity of the gas. We 
plan to study this elsewhere. In addition, our study in- 
dicates that the highly magnetized gas in gravitational 
potential is subject to instabilities (Parker-Rayleigh- 
Taylor-type) which drive turbulence and induce recon- 
nection diffusion of magnetic field. This is another av- 
enue that we intend to explore. 

We report fast magnetic diffusion which happens at 
the rate of turbulent diffusion, but within the present 
set of simulations we do not attempt to precisely evalu- 
ate the rate. Thus we do not attempt to test, e.g., the 
predictions in iLazarianl (|2006| ) of the variations of the 
turbulent diffusion rate with the fluid magnetization for 
the passive scalar field. We also observe that while the 
magnetic field and the passive scalar field diffuse fast, 
there are differences in their diffusion arising, e.g., from 
magnetic field being associated with magnetic pressure. 
We have not attempted to quantify these differences in 
our work either. 

The justification of our results being applicable to 
molecular clouds is based on the model of fast mag- 
netic reconnection in the partially ionized gas in 
lLazarian et al.l (|2004fi . This model is a natural general- 
ization of the LV99 model, but, unlike the LV99 model, 
it has not been numerically tested yet. Such testing will 
be valuable in view of our present study. 



9. SUMMARY 

Motivated by a vital problem of the dynamics of mag- 
netic fields in astrophysical fluids, in particular, by the 
magnetic flux removal in star formation, in this paper we 
have numerically studied the diffusion of magnetic field 
both in the absence and in the presence of gravitational 
potential. Recent work on validating the idea of LV99 
model of reconnection supports our assertion that our 
results obtained at moderate resolution represent the dy- 
namics of turbulent magnetic field lines in astrophysics. 
Our findings obtained on the basis of three-dimensional 
MHD numerical simulations can be briefly summarized 
as follows: 

1. In the absence of gravitational potential the "re- 
connection diffusion" removes strong anti-correlations of 
magnetic field and density that we impose at the start of 
our simulations. The system after several turbulent eddy 
turnover times relaxes to a state with no clear correla- 
tion between magnetic field and densi ty, reminiscent of 
the ob servations of the diffuse ISM by iTroland fc Heilesl 

(psl . 

2. Our simulations that started with a quasi-static 
equilibrium in the presence of a gravitational potential, 
revealed that the turbulent diffusivity induces gas to con- 
centrate at the center of the gravitational potential, while 
the magnetic field is efficiently pushed to the periphery. 
Thus the effect of the magnetic flux removal from collaps- 
ing clouds and cores, which is usually attributed to am- 
bipolar diffusion effect, can be successfully accomplished 
without ambipolar diffusion, but in the presence of tur- 
bulence. 

3. Our simulations that started in a state of dynamical 
collapse induced by an external gravitational potential 
showed that in the absence of turbulence, the flux-to- 
mass ratio is preserved for the collapsing gas. On the 
other hand, in the presence of turbulence, fast removal 
of magnetic field from the center of the gravitational po- 
tential occurs. This may explain the low magnetic flux- 
to-mass ratio observed in stars compared to the corre- 
sponding ratio of the interstellar gas. 

4. As an enhanced Ohmic resistivity to remove mag- 
netic flux from cores and accreti on disks has been ap- 
pealed in the literature, e.g., by IShu et ail (|2006l ). we 
have also compared models with a turbulent fluid and 
models without turbulence but with substantially en- 
hanced Ohmic diffusivity. We have shown that, in terms 
of the magnetic flux removal, the reconnection diffusion 
can mimic the effect of an enhanced Ohmic resistivity. 

5. In addition, our results extend earlier findings of 
ICho et all ()2003al ) for heat advection by magnetized tur- 
bulence. We show that heat advection can be parameter- 
ized by the product of the turbulence injection scale and 
the turbulent velocity for a range of Alfvenic and sonic 
Mach numbers. 
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APPENDIX 



A. TEST CASE: TURBULENT ADVECTION OF HEAT 

We present here some studies of three-dimensional 
transport properties of M HD turbulence. O ur setup is 
identical to that adopted in lCho et al.l ((200 3a). Our goal 
is to provide simulations with higher numerical resolution 
(the resolution employed in iCho et all 12003a! is mostly 
192 3 ) and for a larger parameter space. In the present 
work, we study how compressibility and magnetization 
affect the results. We note that at the time when the 
(|2005ah 



results in lCho et al.1 ()2003aj ) were obtained the issues of 
reconnection in turbulent media were more speculative. 
Thus, after nu merical testing of th e model of reconnec- 
tion in LV99 by iKowal etalT p009) . we feel it is appro- 
priate to revisit that domain. 

Besides the ideal MHD set of equations, we also evolve 
four independent fields of passive scalars $ct(x) (a — 
1, 2, 3, 4) by the continuity equation: 



dt 



■V-(* Q u) = 0. 



(Al) 



These passive scalar fields can trace, for example, the vol- 
ume concentration of some physical property of the gas, 
like metalicity or heat 12 , as long as the time-scale associ- 
ated with the molecular (microscopic) diffusion is larger 
than the typical dynamical time-scale of the flow. In our 
simulations this is ensured by fluids being turbulent 13 . 

These passive scalar fields are assumed to have initial 
spherical symmetry with a Gaussian radial profile: 



* a (x,t = 0) 



expi- 



3(x-x ) 2 



(A2) 



where x is the center of the box and, as in ICho et all 
(|2003aD . we choose the initial dispersion er = Ly/3/l6V% 
(« 19 grid cells, for the resolution employed of 256 3 ), to 
ensure that the characteristic width ero is inside the in- 
ertial range of the turbulence. This initial distribution 
of <fr a is a natural choice to study its diffusion. If we 
had, for instance, a microscopic diffusion coefficient k, 
the Gaussian shape of the distribution would remain un- 
altered, with the dispersion increasing linearly in time at 
a rate k. 

The initial magnetic field Bo is uniform and parallel 
to the x-direction, and the initial density field is also 
uniform with p = 1. The passive scalar fields are injected 
after the turbulence is fully developed. The first scalar- 
field $i is injected at t = 3 (approximately 6 turn-over 
times), and the other fields (a = 2,3,4) at each time step 
after. 

12 It may appear that there is some contradiction here between 
the assumption of an isothermal equation of state (which implies 
instantaneous heat diffusion) and the adoption of a variable scalar 
field to describe heat. However, if we had employed an adiabatic 
equation of state, the results below would not change se nsitively, as 
long as turbulent transport of heating is concerned fe.e.. lCho et al.l 
2003a). 

13 As soon as the molecular diffusivity rate L 2 /rj gets larger than 
the rates associated with the flow, i.e., V/L, the corresponding 
Reynolds number Re = VL/r] gets less than unity and turbulence 
decays at the scale L. 



Table [5] shows the set of parameters studied. The mag- 
netic field is given in terms of the initial Alfven speed. 
Turbulence is injected with the same power e = 1 in 
all the models. These runs cover four combinations of 
regimes of sonic and Alfvenic Mach numbers (M s and 
Ma, respectively). The quantities shown in Table[5](Afg, 
Ma, Vi- ms , -Brms) are averages taken over the time after 
the injection of the first passive field (these time averages 
are over the available computed cubic domains taken at 
every 0.25 time step). The standard deviation of these 
averages are shown within parentheses. For each data 
cube, the Mach numbers correspond to the average of 
the local values computed over the entire box. The res- 
olution employed in all the simulations is 256 3 . 

A.l. Results 

Compressibility can change the properties of flows sub- 
stantially. For instance, high Mach number turbulence is 
known to create regions of enhanced density which will 
coexist with the expanses of low density. At the same 
time, one should expect to see similarity between the 
properties of fluid in low-Mach number flows and incom- 
pressible flows. 

We first consider a case of heat transport, which can be 
considered as a test case, as we can compare our results 
with th e earlier heat transport simulations in Ch o et al.l 
(|2003al) . This is the case where the back-reaction can 
be ignored and we can use a passive scalar diffusion to 
represent the diffusion of heat. 

Let us relate the evolution of the dispersion a with 
the turbulent diffusion coefficient ryturb, that is, the co- 
efficient that gives the rate at which the passive scalar 
field diffuses in scales larger than the turbulent scale. 
The dispersion a is calculated through the definition: 



where x 



/ x<£(x)d 3 x 



If A is the 



rms distance between two fluid elements being advected, 
for A higher than the injection scale of the turbulence /i n j , 
its evolution is related to ?7turb by 5 A 2 ~ rjturbSt. 

Considering ordinary hydrodynamic turbulence, we 
can write, for the rms distance I between two fluid ele- 
ments, within the inertial range, 



61 ~ vilSt, 



(A3) 



where vi is the velocity at scale I. From the assumption 
that both relations above should be valid at the scale /i n j , 
it must be true that ryturb ~ ^injWturb, and we could define 
a constant Cd yn so that: 

f?turb = Cdyn'inj^turb- (A4) 

Now, one wonder if Cd yn is reduced by the magnetic 
field in magnetized turbulence (due to suppression of tur- 
bulent mixture). Cdyn is rela ted to the constant of pro- 
portionality in equation (|A3j) . Again, let I be the dis- 
tance between two fluid elements initially separated by 
a distance within the inertial range in ordinary hydro- 
dynamic turbulence. Repeating equation (|A3[) and using 
the Kolmogorov's phenomenology, 



51 ~ viSt ~ {el) 1/3 St, 



(A5) 
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TABLE 5 

Simulations of the Turbulent Diffusion of Passive Scalar Fields 



Model 


c s 


Bo 


M s 


M A 


Kms 


^rms 


c» 


c»/v ms 


Resolution 


Al 


3.0 


0.1 


0.25(0.00) 


2.41(0.06) 


0.75(0.01) 


0.51(0.01) 


0.4 


ps 0.5 


256 3 


A2 


3.0 


1.0 


0.27(0.01) 


0.78(0.03) 


0.80(0.02) 


1.23(0.01) 


0.15 


ps 0.2 


256 3 


A3 


0.1 


0.1 


7.53(0.13) 


2.13(0.12) 


0.75(0.01) 


0.38(0.01) 


0.4 


ps 0.5 


256 3 


Al 


0.1 


1.0 


7.44(0.12) 


0.55(0.01) 


0.74(0.01) 


1.14(0.01) 


0.3 


w 0.4 


256 3 



where e is the power of injection (or dissipation) of the 
kinetic energy. Integrating the last expression, 

lV3-f /3 = C R eV 3 (t-t ), (A6) 

where the dimensionless constant Cr is th e Richardson 
constant (Richardson 1926; iLesieurl [1990D . Therefore, 
Cdyn is related to Cr. 

Using the dispersion a of our experiments as the rms 
distance I in the equation (A6), and supposing that this 
equation remains valid in MHD turbulence, 

ff 2/3_ a 2/3 = cu l/3 (t _ to)) (A7) 

where C* has not necessarily the same value as that of 
Cr. 

Figure H31 shows the evolution of (it 2 / 3 — cr^/ 3 ) for our 
models. We have made an offset of the values for each 
model in order to make the visualization easier. Table 
[5] lists the approximated values of C* for the models, 
obtained from a linear fit of the data. 



1.6 
1.4 


A1: 
A2: 
A3: 
A4: 


M s ~0.2, M A ~2.4 

Mc-0.3, M A ~0.8 

M s =7.5, M A =2.1 

M s =7.4, M A ~0.5 — — 


1.2 






1.0 






0.8 






0.6 




j*^;***^ 
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Fig. 13. — Evolution of <r 2 / 3 — a^J 3 for different models, where 

2 f (x-x) 2 $(x)d 3 a; , _ fx$(x)d 3 i , . , , , 

a = ^ — r , ' .).,' and x = ■' r ' s — . We nave introduced an 

j ^{xjd^x j "^(xjd^a; 

offset for the lines of each model in order to make the visualization 
clearer. 

All in all, our results presented in Figure [TS] arc con- 
sistent with the scalings of diffusion shown in Ch o et all 
(|2003aft . This indicates that the relation (|A7|) contin- 
ues approximately valid in the MHD case. The slopes 



do not seem to be very sensitive to M s , but we can 
note a small difference for the trans- Alfvenic and super- 
Alfvenic cases, the constant C* being a little higher in 
the super- Alfvenic cases. This finding is consistent with 
the prediction of changes of the diffusion efficiencies of 
ILazarianl (|2006D . There it is shown that for higher mag- 
netization, corresponding to an Alfven Mach number 
Ma = V/Va < 1, where V is the injection velocity and 
Va is the Alfvenic velocity, the turbulence at the large 
scales is weak, which produces weak mixing. In our pa- 
rameterization this corresponds to the decrease of C*, 
which is the case of A2 simulation. Interestingly enough, 
in the case of Ma < 1 and M s ^> 1, which is the case 
of A4 simulation, we see the increase of C*. T he case of 
high M ach number turbulence is not covered by ILazarianl 
(2006), but we expect that the concentration of gas in 
dense regions plays an important role for the effect. In 
regions of gas concentration the Alfveni c velocity drops 
and t he corresponding Ma increases (seefBurkhar t et al.l 
2009). As a result, the mixing should be less constrained 
by the magnetic field as we see in the A4 model. Compar- 
ing cases Al and A3 we observe that for super- Alfvenic 
turbulence and different Mach number M s the diffusivity 
is rather sim ilar. 

Note that iCho et all (|2003al ) obtained the values of 
Ct/Vrns = 0.4 (V A = 0.1 and M s = 2.3), 0.3 (V A = 1 
and M s = 2.3), 0.4 (V A = 1 and M s = 0.3), and 0.5 
(hydrodynamic case and M s = 0.3), wh ich are consistent 
with our findings. In lCho et al1(|2003al) . V rms » 0.8 in all 
the cases (to make a comparison we had to convert the 
time-scale to the present code units, i.e. 1 time step here 
= 2-7T time steps there). We conclude that, for a wide 
range of magnetizations and Mach numbers, turbulence 
is efficient in inducing turbulent diffusivity of heat. 

Our results show that the turbulent diffusivity of pas- 
sive scalar fields r?turb is well described by ryturb = 
C*^turb^inj, similar to the hydrodynamic case, and that 
the coefficient C* is not very sensitive to the strength 
of the magnetic field. Therefore, the magnetic field does 
not impose a strong suppression of the turbulent diffusion 
perpendicular to it as far as heat transfer is c oncerned. 

As already pointed in ICho et all (|2003af ) (see also 
ILazarianl I2006T ). the effectiveness of the turbulent diffu- 
sion has important consequences on the thermal diffusion 
in the ISM and intracluster medium (ICM). The nearly 
isotropic turbulent diffusivity can be 1 order of magni- 
tude higher in the gas of central regions of galaxy clus- 
ters, like the Hydra A cluster, compared to the Spitzer 
value, based on the kinetic theory. In the ISM, in mix- 
ing layers, where turbulence is very strong, the turbulent 
diffusivity can be 2 orders of magnitude higher than the 
laminar values (see lEsauivel et aLll2006|) . Our results ex- 
tend the range of the applicability of the turbulent heat 
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advection model of iCho et al.l (|2003aft . 

B. DIFFUSION OF PASSIVE SCALAR IN THE ABSENCE 
OF GRAVITY 

In our simulations described in Section 4 (see Table [TJ 
we also keep track of passive scalar field representative 
either of heat or passive impurity. 

Figure [14] (left panels) shows the time evolution of 

(^>0.25/(P)o.25-^/P and ( $ )p.25/(P)o.25-*/P- B ° th 

quantities have a similar behavior, and all models seem 
to achieve the characteristic average values (B z / p and 
$/p) roughly at the same time. 

Right plot of Figure [TJ] is an analog of Figure [3] A com- 
parison between these figures reinforces that the quantity 
(B z ) z I (p) z evolves similar to (<£>) z / (p) z - How can we 
understand this behavior? 

Using the expression ^=Vx(vxB) one gets ^ = 
(B-V)v— (v-V)B, when V-B = is accounted for. Using 
the continuity equation in the form V • v 
one can get 



l dp 



p dt 



which coincides with Equation lAll if $ Q is substituted by 
B/p. 

Figure [15] (left panels) shows the distribution of (pK 
versus ($) z for the model B2 (Bq = 1.0, see Table \U 
at different time steps (t ~ 0, 10). The initial anti- 
correlation gives place to a tight correlation between the 
passive scalar and density. Right side of Figure [15] shows 
the evolution of (<5<I>, 8p) (defined by equation where 
B must be substituited by $). 

C. FURTHER COMPARISON OF MAGNETIC 
DIFFUSIVITY AND RESISTIVITY EFFECTS IN THE 
PRESENCE OF GRAVITY 

Figure \W\ shows the comparison of the evolution of 
(B Z ) R=0 ? - L for models CI, C3, C4, C5, C6, and C7 
of Table [2] with similar resistive models without turbu- 
lence of Tabled] with different values of Ohmic diffusivity 
VOhm- From these comparisons (and others where (B z 
is measured for different values of R, like in Figure [~~ 
we have estimated the values of ?yturb shown in Table [2] 



d „\ B B „ 

— + v • V - = - • Vv 

dt ) p p 



(Bl) 

REFERENCES 



Beresnyak, A., & Lazarian, A. 2006, ApJ, 640, L175 

Beresnyak, A., & Lazarian, A. 2009a, ApJ, 702, 1190 

Beresnyak, A., & Lazarian, A. 2009b, ApJ, 702, 460 

Biskamp, D. 1986, Physics of Fluids, 29, 1520 

Boldyrev, S. 2005, ApJ, 626, L37 

Boldyrev, S. 2006, Phys. Rev. Lett., 96, 115002 

Burkhart, B., Falceta-Gongalves, D., Kowal, G., & Lazarian, A. 

2009, ApJ, 693, 250 
Cho, J., & Lazarian, A. 2002, Phys. Rev. Lett., 88, 245001 
Cho, J., & Lazarian, A. 2003, MNRAS, 345, 325 
Cho, J., Lazarian, A., Honein, A., Knaepen, B., Kassinos, S., & 

Moin, P. 2003a, ApJ, 589, L77 
Cho, J., Lazarian, A., & Vishniac, E. T. 2003b, ApJ, 595, 812 
Dobrowolny, M., Mangeney, A., & Veltri, P. 1980, Phys. Rev. 

Lett., 45, 144 
Draine, B. T., & Lazarian, A. 1998, ApJ, 508, 157 
Elmegreen, B. G., & Scalo, J. 2004, ARA&A, 42, 211 
Esquivel, A., Benjamin, R. A., Lazarian, A., Cho, J., & Leitner, 

S. N. 2006, ApJ, 648, 1043 
Falceta-Goncalves, D., Lazarian, A., & Kowal, G. 2008, ApJ, 679, 

537 

Galli, D., Lizano, S., Shu, F. H., & Allen, A. 2006, ApJ, 647, 374 

Giacalone, J., & Jokipii, J. R. 1999, ApJ, 520, 204 

Gogoberidze, G. 2007, Phys. Plasmas, 14, 022304 

Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763 

Heitsch, F., Zweibel, E. G., Slyz, A. D., & Devriendt, J. E. G. 

2004, ApJ, 603, 165 
Higdon, J. C. 1984, ApJ, 285, 109 
Iroshnikov, P. S. 1963, AZh, 40, 742 

Kowal, G., Lazarian, A., & Beresnyak, A. 2007, ApJ, 658, 423 
Kowal, G., Lazarian, A., Vishniac, E. T., & Otmianowska-Mazur, 

K. 2009, ApJ, 700, 63 
Kraichnan, R. H. 1965, Phys. Fluids, 8, 1385 

Kulpa-Dybel, K., Kowal, G., Otmianowska-Mazur, K., Lazarian, 
A., & Vishniac, E. 2009, A& A, in press l larXiv:0909.1265l l 

Lazarian, A. 2005, in AIP Conf. Proc. 784, Magnetic Fields in the 
Universe: From Laboratory and Stars to Primordial Structures 
(Melville, NY: AIP), 42 

Lazarian, A. 2006, ApJ, 645, L25 



Lazarian, A., & Beresnyak, A. 2006, MNRAS, 373, 1195 
Lazarian, A., & Vishniac, E. T. 1999, ApJ, 517, 700 
Lazarian, A., & Vishniac, E. T. 2009, RevMexAA, 36, 81 
Lazarian, A., Vishniac, E. T., & Cho, J. 2004, ApJ, 603, 180 
Leao, M. R. M., de Gouveia Dal Pino, E. M., Falceta-Goncalves, 

D., Melioli, C., & Geraissate, F. G. 2009, MNRAS, 394, 157 
Lesieur, M. 1990, NASA STI/Recon Technical Report A, 91, 

24106 

McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565 
Mestel, L., & Paris, R. B. 1984, A&A, 136, 98 
Mestel, L., & Ray, T. P. 1985, MNRAS, 212, 275 
Mestel, L., & Spitzer, L., Jr. 1956, MNRAS, 116, 503 
Montgomery, D., & Turner, L. 1981, Physics of Fluids, 24, 825 
Padoan, P., Jimenez, R., Juvela, M., & Nordlund, A. 2004, ApJ, 

604, L49 
Parker, E. N. 1958, ApJ, 128, 664 

Passot, T., & Vazquez-Semadeni, E. 2003, A&A, 398, 845 
Petschek, H. E. 1964, In Procs. AAS-NASA Symposium, The 

Physics of Solar Flares, Eds. W. N. Hess. Washington, NASA 

Special Publication, 1964, 50, 425 
Richardson, L. F. 1926, Proc. R. Soc. A, 110, 709 
Shay, M. A., Drake, J. F., Denton, R. E., & Biskamp, D. 1998, 

J. Geophys. Res., 103, 9165 
Shay, M. A., Drake, J. F., Swisdak, M., & Rogers, B. N. 2004, 

Phys. of Plasmas, 11, 2199 
Shebalin, J. V., Matthaeus, W. H., k, Montgomery, D. 1983, J. 

Plasma Phys., 29, 525 
Shu, F. H. 1983, ApJ, 273, 202 

Shu, F. H., Galli, D., Lizano, S., & Cai, M. 2006, ApJ, 647, 382 
Sweet, P. A. 1958, Electromagn. Phenom. Cosm. Phys., 6, 123 
Tassis, K., & Mouschovias, T. C. 2005, ApJ, 618, 769 
Troland, T. H., & Heiles, C. 1986, ApJ, 301, 339 
Uzdensky, D. A., & Kulsrud, R. M. 2000, Phys. of Plasmas, 7, 
4018 

Yamada, M., Ren, Y., Ji, H., Breslau, J., Gerhardt, S., Kulsrud, 

R., & Kuritsyn, A. 2006, Phys. Plasmas, 13, 052119 
Zweibel, E. G. 2002, ApJ, 567, 962 



22 



Santos-Lima et al. 




time time time 

Fig. 14. — Left: evolution of the ratio of the averaged magnetic field over the averaged density (more left) and of the ratio of the averaged 
passive scalar over the averaged density (more right) within a distance R = 0.25L from the central z-axis. The values have been subtracted 
from their characteristic values B/p in the box. Right: evolution of the rms amplitude of the Fourier modes (k x , k y ) = (±1, ±1) of ($) 
(upper curves) and (3>) / (p) (lower curves). The curves for <P were multiplied by a factor of 10. All the curves were smoothed to make 
the visualization clearer. 
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Fig. 15. — Left: distribution of (p) vs. (3>) z for model B2 (see Table[T}, at t = (most left) and t = 10 (most right). Right: correlation 
between fluctuations of the passive scalar field ((5 < 1?) and density (<5p). 




